Aperiodic Quantum Random Walks 
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We generalize the quantum random walk (QRW) protocol for a particle in a one-dimensional chain, 
by using several types of biased quantum coins, arranged in aperiodic sequences, in a manner that 
leads to a rich variety of possible wave function evolutions. Quasiperiodic sequences, following the 
Fibonacci prescription, are of particular interest, leading to a sub-ballistic wavefunction spreading. 
In contrast, random sequences leads to diffusive spreading, similar to the classical random walk 
behaviour. We also describe how to experimentally implement these aperiodic sequences. 

PACS numbers: 



A quantum random walk (QRW) is a natural extension 
to the quantum world of the ubiquitous classical random 
walk. It was first proposed in |lj and widely investigated 
recently (see the recent review by Kempe |2j), mostly in 
connection with possible applications to quantum algo- 
rithms |3|, y|. The generic discrete QRW consists in a 
particle moving on a graph , in a direction depending 
on its internal state (either called spin or chirality), the 
simplest case being a spin i particle on a periodic chain. 
In between each moving step, an unitary transformation, 
called a "quantum coin operator" (QCO), is acted on the 
particle spin state and shuffles the spin related amplitude 
of the particle wave function. Most studies have been fo- 
cused on the so-called Hadamard transform, but more 
general QCO can been used. A main difference between 
the classical and quantum walks is seen on the particle 
spreading, as measured by the long time dependence of 

the standard deviation a(t) = J [x 2 ) t — (x) t . The clas- 
sical case displays a diffusive behavior ( a(t) ~ £ = ) while 
the quantum case is ballistic (er(£) ~ t), as can be proved 
in one dimension from the exactly computed solution [£| . 
The latter result relies on the space periodicity of the 
QRW process, which allows for a Fourier transformed 
wave function simple form. This reminds the behavior of 
tight-binding Bloch electrons under standard quantum 
evolution on a periodic lattice. It is therefore tempting 
to check whether well-known effects of quasiperiodicity 
in the latter case (like a sub-ballistic scaling with time of 
the standard deviation|6j, or the auto-correlation func- 
tion J3) can also be observed in QRW. We address this 
question here by generalising the QRW to the case where 
different quantum coins are applied along three types of 
sequences, either periodic, quasi periodic and random. 



The particle displacement is along a one dimensional 
periodic chain indexed by k € Z, with a corresponding 
orthonormal basis {|fc)} spanning the position Hilbcrt 
space Hp,. To the quantum coin part corresponds a two- 
dimensional Hilbert space Hq spanned by {||) , ||)}- The 



particle wave function reads 

|*)=53o(fl,fc)|*)®|fc) with sG {T,|} (1) 

s,k 

The QRW unitary step operator S(a) is the concatena- 
tion of a displacement D which reads 

D = J2 (IT) (Tl ® \k + 1} (k\ + |4) (II ® I* - 1) (k\) (2) 

k 

following a quantum coin operator C(a), a unitary 
transformation acting on the spin sector. Here we mainly 
use the rather simple QCO, 



C(a) 



cos(a) sin(a) 
sin(a) — cos(a) 



(3) 



which depends on a single parameter a £ [0, ^] . Note 
that C(a) squares to the identity transformation. 

For a — 7r/4, one recovers the widely studied 
Hadamard walk, with its ballistic behavior. Other val- 
ues of a also also gives rise to the same behavior, with 
different prefactors, except for a = ir/2 (corresponding 
to the <t x Pauli matrix) for which the particle remains 
confined. Note that for a = 0, the QCO reduces to the 
diagonal Pauli matrix a z . In this case, the motion is 
completly decoupled into ballistic (right going) spin up 
and (left going) spin down parts, the latter acquiring a 
■k phase at each step. Figure ^ displays the typical wave 
function speading for several a values, with an initial ket 
(IT) +i Il))/v2) which allows for a symmetric probability 
distribution. 

Let us now combine two different step operators S(a) 
and S(f3) into infinite sequences and compute the wave 
function spreading as measured by the exponent c(a, (3) 
in a ~ ^(".Z 3 ). For periodic sequences, the long time 
behaviour is still found to be ballistic (see below for 
quasiperiodic approximants), although displaying a more 
complex structure at the scale of one period. This sug- 
gests that new behaviours can be expected when the pe- 
riod size tends to infinity, as for the quasiperiodic case 
that we now study. 
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FIG. 1: Probability distributions after 100 time steps for 10 
equally spaced a values between and 7r/2. The middle curve 
corresponds to the standard Hadamard case. For a — 0, the 
QCO is simply a z , and the quantum walk amounts to two 
unrelated, left-right, ballistic moves, whose probability distri- 
bution is not distinguished from the bounding box vertical 
edges. For a = n/2, the QCO is a x , and the quantum walk 
is confined near the origin 



We consider a Fibonacci sequence, obtained by iter- 
ation of the recursive rule S n +i = S n S n ~i, with So — 
S(a) and Si = S(a)S(f3). Given a sequence S n , the 
next one can also be otained using a substitution rule 
S(a) -> S(a)S(/3), S((3) -> S(a). Only writing the a and 
(3 symbols, the first sequences, called approximant se- 
quences, read: a, a/3, a(3a, a(3aaf3, af3aal3a(3a . . .. The 
infinite sequence is not periodic. Indeed the occurence 
ratio of a versus tends to the irrational golden mean 
t = (1 + \/5)/2, which cannot be displayed in the re- 
peated unit cell of a periodic sequence. This sequence 
is well known to be quasiperiodic, a kind of order that 
has been widely investigated in the last 20 years in the 
context of quasicrystals. The effect of a sequence of 
quasiperiodic unitary tranformations applied to a spin 
1/2 has been studied by Sutherland |8| and displayed a 
rich behaviour. The latter transformations were more 
generic than the simple QCO used here, but were not 
coupled to a displacement, as in the QRW. In the space 
sector, this sequence, when coding a quasiperioc poten- 
tial, is well known to cause a sub-ballistic behaviour for 
tight-binding electrons |(|,|Z|. In contrast, periodic rep- 
etition of the approximant sequences eventually leads to 
a ballistic spreading. 

A similar behaviour is found here for the quasiperiod- 
ically shuffled QRW (with generic values of a and (3). 
Recall that the position space is here periodic, and that 
the quasiperiodic modulation is applied with time in the 




FIG. 2: Standard deviation a(t) of the probability distribu- 
tions for periodically repeated approximant sequences (with 
period length of respectively 2,3,5 and 8), and for the asymp- 
totic Fibonacci quasiperiodic sequence (grey colour). The lat- 
ter clearly displays a sub-ballistic slope. These curves are 
obtained with the parameters (a = 7r/3, /3 = 7r/6) 



spin sector. The simulation was done for typically sev- 
eral thousands random walk steps. Let us first com- 
pare the spreading with time for periodically repeated 
small approximants and for the asymptotic quasiperi- 
odic sequence. By the latter, we mean that we numer- 
ically generate a long sequence whose length is larger 
than both the chain length and the number of random 
walk steps. The standard deviation is plotted in figure 
13 which clearly displays a qualitative difference between 
the (asymptotically linear with time) periodic case and 
the slower quasiperiodic case. 

The sub-ballistic behaviour is generic in the quasiperi- 
odic case, whatever a and (3 are. But the asymptotic 
slope c(a, (3) is not a smooth function, as seen on figure 
13 The "diagonal " a = (3 corresponds to the periodic 
case and therefore to the expected ballistic slope c = 1. 
More surprising, and not yet fully understood, are the 
clearly visible transverse crest lines, whose (equal) inverse 
slope is very close to the golden mean. To check whether 
this apparent arithmetical relation was not a simple con- 
cidence, we tried an alternative quasiperiodic sequence, 
based on the 'silver' mean : 1 + \[2. Parallel crest lines 
whose slope are simply related to the silver mean were 
again found. We are then tempted to appeal to the 
subtle properties of some iterated maps based on these 
quasiperiodic sequences 9] , displaying either periodic (in 
the order of the Fibonacci sequence) or chaotic behaviour 
which are clearly seen in the computed quasiperiodically 
rotated spin system [8J. The above Fibonacci quantum 
coin sequence proves to have a simple cyclic behaviour 
for any value of the pair (a, (3). But these crest lines also 
appear if we replace the aboce QCO (expression Q) by a 
simple 2 dimensional rotation matrix, in which case the 
quasiperiodically rotated spin system is not generically 
cyclic. 



Note that the QRW possible experimental impleman- 
tation (see below) is of high interest in the context 
of quasiperiodic systems, since, in that case, the sub- 
ballistic behaviour, although clear from many computa- 
tions, has not been yet clearly demonstrated in experi- 
ments. 
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FIG. 3: Slope of the standard deviation a(t), versus the pa- 
rameters (a, j3) for a Fibonacci sequence. 

Let us finally consider random sequences. We are still 
willing to compare with tight-binding electron evolution. 
In the latter case, a one-dimensional random potential 
is expected to generate a (very) long time wave packet 
localisation, preceeded, at short time, by a ballistic mo- 
tion whose range depends on the error range width. We 
consider here two types of disordered QRW. We first gen- 
erate 50 — 50 random sequences with two different QCO, 
defined by a, n/2 — a. To each random sequence corre- 
sponds a definite wave-function spreading, and we there- 
fore average over many disorder realisations to check for 
an asymptotic regime. The second case under consider- 
ation is that of a continuous set of QCO, whose distri- 
bution, centred on a — vr/4, has variable width. This 
situation is related, but not identical, to a model of er- 
rors for an experimental implementation, which has been 
considered in [10( for periodic QRW in an optical lattice. 
In both cases, a diffuse regime, c ~ 0.5 is found, and the 
probability distribution have a gaussian-like shape with 
short range structure (figure 0J. 

Experimental implementation: We can now discuss 
how the above computed behaviours could be experi- 
mentally testet. There are, up to now, several experi- 
mental proposals for the realization of the usual QRW, 
among which, for example, one atom in an optical lattice 
[l0| . trapped ions [l^. systems usinglinear optics [b| 
and cavity quantum electrodynamics [jj, [ll( . A continu- 
ous time version of the algorithm performed in a circular 
chain has already been realized using nuclear magnetic 
resonance 15] and the difference with the classical ran- 




FIG. 4: Typical probability distribution in the '"continuous"' 
random case. A " 'hairy" '-lige gaussian shape is found, with 
a diffusive (c ^ 0.5) standard deviation 



dom walk demonstrated. The realization of the aperiodic 
QRW involving biased alternating coins, represent only 
limited modifications to most of these protocols, imply- 
ing that there is no major difficulty to the experimen- 
tal implementation of the present aperiodic two coins se- 
quences. We can study the specific example of optical 
lattices and describe how the protocol of aperiodic QRW 
can be implemented therein. In the proposal presented 
in [lfj , the internal states of a Rb neutral atom (consid- 
ered as a two state system) are subject to the quantum 
coin. This atom is located initially at one site of an one 
dimensional lattice, and will move to one of its neighbour- 
ing site according to its internal state. The conditional 
translation is performed in the same way as proposed 
in 14] : different internal states feel different optical po- 
tentials and they are kept in the ground state of their 
respective potential. For the case studied in [ljj, the 
relevant internal atomic states are the hyperfine struc- 
ture states \F = l,m/ = 1) and \F — 2, m/ = 2). They 
will be denoted from now on as qubits |0) and |1). The 
corresponding potential is Vo(x,8) = [V mg= i/ 2 (x, 8) + 
3V m .=_i/ 2 (M)]/4 and V x {x,9) = V m ^ /2 (x,9), where 



V m ,=±i/ 2 (x,6) = a\E \ 2 sin (kx± 6) |14|. The angle 9 
is half the angle between the polarizations of the two 
lasers that form the lattice and E is the amplitude of 
the electrical field. For the specific choice of 9 = 0, the 
minima for both hyperfine states coincide. The parame- 
ter 9 can be adiabatically changed by turning one of the 
lasers polarization. The potential wells corresponding to 
each one of the internal states are thus translated with 
respect to each other. Since the minima for both inter- 
nal states coincide again at 9 = n, at this point the two 
possible internal states are present in the same position. 
The "mixing" step, corresponding to the Hadamard gate, 
is performed by a laser tunned to the frequency separat- 
ing the two atomic internal states. The application of 
the laser pulse occurs when the minima for both states 
coincide, and its duration determines the superposition 



of the two internal states that is created. The hamilto- 
nian corresponding to the atom-laser interaction can be 
written as: 



F int = n(|0)(l|e^ + |l)(0|e-^, 



(4) 



where fi is the Rabi frequency, e 1 ^ is the pulse's phase. 
The time evolution reads : 



i(j> 



|0> {t) = cos (fit) |0) - sin (fit)e 
|1> (t) = sin (fit)e^ |0) + cos (fit) 



1) 



ID 



(5) 



This means that an appropriate choice of the laser's phase 
and of the pulse duration can build all possible superpo- 
sitions of the two states. 

The modifications to the usual QRW protocol sug- 
gested in this paper demand changing only this point 
of the experimental proposal: the generalized coins cre- 
ate biased superpositions of the internal states, and their 
realization is possible by controlling the laser pulse dura- 
tion. In order to alternate two different GCO, one only 
needs to alternate two different pulse durations. 

Conclusion: We have shown how the standard Quan- 
tum Random Walk framework can be considerably en- 
riched by allowing more than one type of Quantum Coin 
Operator, arranged along different squences. In particu- 
lar, quasiperiodic binary sequences lead to a sub-ballistic 
wave packet spreading characterised by sequence depen- 
dant slopes. The quantum state evolution depends on 
both the precise values taken by the two coins and the 
sequence itself. Note that we have studied here the sim- 
plest cases, with Fibonacci sequences and generalized 
Hadamard coins. More complex evolution are expected if 
the two coins are picked in a more generic set, or arranged 
along more complex (even not random) sequences. We 
have also studied random sequences, which leads to dif- 
fuse spreading, in contrast with the localisation effect en- 
countered for quantum particles subjected to disordered 
potentials in one dimension. 

A very interesting aspect of these binary quantum ran- 
dom walk is their possible experimental implementation. 



Indeed sub-balistic spreadings are often computed (spe- 
cially in the quasiperiodic case), but rarely observed. We 
have described here the mofication implied by going from 
single to binary sequences, which is not in principle very 
complicated. 
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